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1.  Introduction 

The  work  presented  in  this  report  is  motivated  by  the  need  to  calibrate  a  new  pul sed-micro wave 
breast  tumor  detection  system  for  patient-specific  skin  parameters.  A  two-dimensional  time- 
domain  inverse-scattering  algorithm  based  upon  the  finite-difference  time-domain  method  is 
presented  for  determining  the  skin  thickness  and  the  relative  permittivity  Er-skm  and  electric 
conductivity  cSki„  of  skin  in  the  microwave  range.  The  algorithm  traces  a  search  trajectory  in  the 
two-dimensional  (£r-skin,<7skin)  parameter  space.  The  minimal  parameter  estimation  error  along  this 
trajectory  yields  a  set  of  approximate  parameter  values.  It  is  shown  that  the  convergence  of  the 
inverse-scattering  technique  depends  on  the  shape  and  the  duration  of  the  illuminating 
electromagnetic  wave  pulse  chosen  for  the  electrical  parameter  reconstruction.  The  time-domain 
nature  of  the  inverse  algorithm  allows  for  limiting  the  region  of  inversion  using  causality.  Thus, 
when  the  parameters  of  the  skin  are  estimated,  the  skin  thickness  can  be  determined  by 
comparing  the  measurement  with  a  simulated  all-skin  response.  Finally,  the  time-domain 
inverse-scattering  algorithm  is  tested  for  robustness  in  the  presence  of  broadband  Gaussian  noise. 
We  note  that  the  causality  of  the  algorithm  can  then  be  exploited  for  the  underlying  breast  tissue 
parameter  recovery  using  the  same  methodology. 

2.  Body  -  Description  of  the  Research  Accomplishments 

Investigator’s  work  and  the  research  accomplishments  supported  by  the  grant  from  May  1999  - 
May  2001  have  been  described  in  the  previous  two  annual  reports  submitted  to  US  Army 
Medical  Research  and  Materiel  Command  Office.  By  the  end  of  August  2001,  the  investigator 
completed  her  Ph.  D.  requirements  and  was  afterwards  therefore  no  longer  supported  by  the 
grant. 

Cumulative  research  accomplishments  supported  by  the  duration  of  the  award  (May  1999 
-  August  2001)  have  been  summarized  in  her  doctoral  thesis  and  in  a  manuscript  prepared  for 
submission  to  a  major  journal  {IEEE  Transactions  on  Antennas  and  Propagation ).  The  paper 
was  reviewed  and  the  author  is  presently  introducing  minor  modifications  to  comply  with  the 
reviewers’  suggestions.  However,  this  manuscript  describes  in  detail  major  findings  of  the  work 
supported  by  the  award  and  is  therefore  appended  to  this  report  [Appendix  A].  The  manuscript 
first  reviews  the  background  literature  on  dielectric  properties  in  the  microwave  range,  skin 
thickness  in  the  area  of  the  human  breast  and  numerical  inverse-scattering  methods.  Next, 
calculations  are  presented  to  illustrate  the  importance  of  knowledge  of  the  correct  skin  thickness 
for  determining  the  time  delays  needed  for  the  image-formation  signal-processing  algorithm  tied 
to  the  microwave  imaging  system.  Then,  results  are  presented  to  show  the  development  of  two- 
dimensional  time-domain  inverse-scattering  algorithm  for  simultaneous  estimation  of  electrical 
permittivity  and  conductivity  of  the  skin  layer.  This  algorithm  locates  a  search  trajectory  in  the 
parameter  space.  Results  are  then  presented  which  show  that  the  convergence  of  the  search 
trajectory  and  the  convergence  error  depend  on  the  shape  and  the  duration  of  the  impinging 
pulse.  The  estimated  skin  electrical  parameters  then  permit  the  skin  thickness  to  be  determined 
by  comparing  the  measurement  with  a  simulated  all-skin  response.  As  the  final  step,  the 
robustness  of  the  algorithm  is  tested  in  the  presence  of  simulated  Gaussian  noise  for  various 
signal-to-noise  ratios.  For  more  details  and  figures  which  illustrate  major  research  findings, 
please  refer  to  Appendix  A. 
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3.  Key  Research  Accomplishments 


The  investigator,  who  has  electrical  engineering  background,  gained  valuable  medical 
insight  into  the  anatomy  of  the  human  breast,  nature  of  the  main  breast  tumor  categories  and 
the  related  symptoms. 

Numerical  models  confirmed  that  it  is  important  to  determine  accurate,  local  patient-specific 
breast  thickness  in  order  to  properly  calibrate  the  microwave  breast  tumor  and  imaging 
system. 

Motivated  by  the  previously  stated  need  for  the  microwave  imager  calibration,  we  have 
developed  a  two-dimensional  (2-D)  inverse-scattering  finite-difference  time-domain  (FDTD) 
algorithm  for  the  recovery  of  the  electrical  parameters  of  the  skin  (permittivity  Er.skin  and 
electric  conductivity  askjn)  in  the  microwave  range.  This  algorithm: 

-  locates  a  search  trajectory  in  the  (Er-skin,<7Skin)  parameter  space; 

-  determines  the  square-normed  error  along  the  search  trajectory; 

-  locates  the  approximate  parameter  values  for  the  minimum  error. 

The  algorithm  was  tested  for  several  excitation  waveforms. 

For  each  waveform  under  investigation,  the  algorithm  was  tested  for  robustness  in  the 
presence  of  zero-mean  additive  Gaussian  noise  for  several  signal-to-noise  ratio  values. 

The  results  imply  that: 

-  the  robustness  of  the  time-domain  inverse-scattering  algorithm  for  recovering  the 
electrical  parameters  of  the  skin  depends  on  the  duration  and  shape  of  the  excitation 
waveform.  Specifically,  of  the  excitations  considered,  a  10-ps  differentiated  Gaussian 
pulse  provides  superior  robustness,  even  relative  to  the  ramp  signal  having  the  same  rise¬ 
time,  if  the  time  window  used  for  the  inverse-scattering  procedure  contained  the  entire 
scattered  waveform. 

-  The  use  of  a  time-linear  excitation  yields  a  linear  search  trajectory  in  the  parameter  space 

(Er-skinjt^skin)' 

-  Proper  combination  of  time-linear  and  short  bipolar  pulse  excitation  data  can  yield  an 
efficient  and  robust  search  strategy. 

By  exploiting  causality,  after  skin  electrical  parameters  are  determined,  skin  thickness  can  be 
estimated  by  comparing  the  measurement  and  the  simulated  all-skin  response. 

The  same  methodology  can  then  be  employed  to  determine  the  electrical  parameters  of  the 
underlying  breast  tissue. 


2 


I 


4.  Reportable  outcomes 

In  this  section,  the  investigator  reports  the  outcomes  which  were  supported  by  the  award  for  the 
entire  award  period  (May  1999  -  August  2001).  The  outcomes  which  were  appended  in  the 
previous  reports  are  only  listed,  and  those  which  date  from  after  the  last  annual  report  are 
included  in  the  appendices  as  noted. 

Doctoral  Thesis 

M.  Popovic,  “Microwave  Imaging  of  Breast  Tissues”,  published  March  13,  2002,  United  States 
Copyright  Office,  Registration  Number  TX  5-498-417,  USA.  (Please  see  Appendix  C  for  the 
front  page.) 

Journal  Papers 

M.  Popovic,  A.  Taflove,  “2-D  FDTD  Inverse-Scattering  Scheme  for  Determination  of  Breast 
Skin  Properties  at  Microwave  Frequencies”,  manuscript  submitted  to  IEEE  Trans.  On  Antennas 
and  Propagation,  August  2001.  Manuscript  is  currently  under  revision  in  order  to  comply  with 
the  reviewers’  comments  prior  to  resubmission.  [Appendix  A] 

Refereed  Conference  Proceedings  Papers 

M.  Popovic  and  A.  Taflove,  “Skin  Thickness  and  Dielectric  Parameter  Evaluation  in  the 
Microwave  Range  Using  an  FDTD  Inverse-Scattering  Technique”,  Proceedings  on  2001  URSI 
International  Symposium  on  Electromagnetic  Theory ,  Victoria,  Canada.  May  13-17,  2001. 

Abstracts  and  Presentations 

M.  Popovic  and  A.  Taflove,  “Nonin vasive  Measurement  of  Breast  Skin  Properties  in  the  Microwave 
Range:  Efficiency  Improvement  of  the  2-D  FDTD  Inverse-Scattering  Scheme”,  The  Tenth  Biennial  IEEE 
Conference  on  Electromagnetic  Field  Computation,  Perugia,  Italy,  June  16-19,  2002.  [Appendix  B] 

M.  Popovic  and  A.  Taflove,  “2-D  FDTD  Inverse-Scattering  Scheme  for  Determining  Microwave  Breast 
Skin  Properties:  the  Results  and  the  Implied  Optimization  Options”,  The  Bioelectromagnetics  Society  24th 
Annual  Meeting,  Quebec  City,  Canada,  June  23-27,  2002.  [Appendix  B] 

M.  Popovic  and  A.  Taflove,  “Skin  Thickness  and  Dielectric  Parameter  Evaluation  in  the  Microwave 
Range  Using  an  FDTD  Inverse-Scattering  Technique”,  2001  URSI  International  Symposium  on 
Electromagnetic  Theory,  Victoria,  Canada.  May  13-17,  2001. 

M.  Popovic  and  A.  Taflove,  “Time  Domain  Inverse-Scattering  Technique  for  Obtaining  Microwave 
Properties  of  Near-Surface  Body  Tissues”,  Progress  in  Electromagnetics  Research  Symposium  (PIERS), 
July  5-14,  2000,  Cambridge,  Massachusetts,  USA,  pp.  550. 

M.  Popovic  and  A.  Taflove,  “Obtaining  Microwave  Properties  of  Near-Surface  Body  Tissues  Using  2-D 
FDTD  Inverse-Scattering  Technique”,  Bioelectromagnetics  Society  (BEMS)  2000  Twenty-Second 
Annual  Meeting,  June  11-16,  2000,  Munich,  Germany,  pp.23. 

S.  C.  Hagness,  M.  Popovic,  A.  Taflove  and  Jack  E.  Bridges,  “An  Image  Reconstruction  Algorithm  for 
Ultra  wideband  Microwave  Radar  Technology  Applied  to  Breast  Cancer  Detection”,  1999  IEEE  AP-S 
International  Symposium  and  USNC/URSI  National  Radio  Science  Meeting,  July  11-16,  1999,  Orlando, 
Florida,  USA. 

Degrees  Obtained 

The  applicant  defended  her  Doctoral  work  in  July  2002,  and  completed  the  requirements  for  the 
Ph.D.  in  August  2002  [Appendix  C],  obtaining  the  official  diploma  in  December  2002. 
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Awards 

-  May  1999:  Award  for  one  of  the  two  best  poster  presentations,  ECE  Research  and  Industry 
Fair,  ECE  Department,  Northwestern  University. 

-  BEMS  2000  Twenty-second  Annual  Meeting,  EMF  Therapeutics  Clinical  Application  Award, 
June  11-16,  Munich,  Germany. 

Employment  and  Research  Opportunities 

The  doctoral  research  funded  by  the  award  was  referred  to  in  over  ten  applications  for  academic 
positions  across  North  America.  One  of  the  interviews  led  to  employment  of  the  investigator  as 
an  assistant  professor  with  the  Department  Electrical  and  Computer  Engineering  at  McGill 
University,  Montreal,  Canada,  beginning  September  2001,  to  the  present  date.  (Please  refer  to 
Appendix  D  for  a  copy  of  the  appointment  letter  from  the  Academic  Vice-Principal.) 

As  an  assistant  professor  and  a  researcher,  the  investigator  applied  for  funding  from  both 
federal  and  provincial  sources  in  Canada  in  order  to  continue  the  research  started  under  the 
Department  of  Defense  pre-doctoral  grant. 

Received  funding: 

•  Natural  Sciences  and  Engineering  Research  Council  of  Canada  (NSERC)  -  Individual 
Research  Grant  (3-year  program;  level  received:  approximately  $42,000  CDN  /  year); 

•  Internal  Research  Grants  Office  Funding  -  McGill  University:  Individual  New  Researcher 
Grant,  $17,500  CDN,  one-time  grant. 

Applied  for  and  pending  funding: 

•  Fonds  pour  la  Formation  de  Chercheurs  et  l’aide  a  la  Recherche  (FCAR)  -  Quebec: 

Programmes  destinees  aux  nouveaux  chercheurs  -  requested  $15,000  CDN; 

-  Programme  strategique  de  professeurs-chercheurs  reserve  aux  universities  -  supplies  the 
salary  of  the  applicant  to  the  university  department  and  reduces  teaching  load  in  favor  of 
the  research  time. 

5.  Conclusions 

In  this  report,  we  presented  the  research  conducted  to  investigate  possibilities  of  non-invasive 
determination  of  near-surface  breast  tissue  electrical  properties  in  order  to  aid  the  patient-specific 
calibration  of  the  microwave  tumor  detection  and  imaging  device.  When  a  layer  of  materials 
(tissues)  is  excited  with  a  signal,  information  about  the  tissue  electrical  properties  is  hidden  in  the 
backscattered  signal  data.  A  2-D  FDTD  inverse-scattering  algorithm  was  developed  to  show  that 
non-invasive  measurement  of  near-surface  tissue  properties  depends  heavily  on  the  shape  and  the 
duration  of  the  impinging  signal.  Moreover,  studies  suggest  that  usage  of  a  short,  bipolar 
excitation  can  yield  an  accurate  estimate  of  the  skin  electrical  properties  even  in  the  presence  of 
high-level  random  Gaussian  noise. 

A  logical  near-term  extension  of  this  work  is  to  construct  an  additional  software-element 
which  would  use  the  same  technique  to  obtain  the  electrical  properties  of  the  breast  tissue 
underlying  the  skin  after  first  determining  the  skin  thickness  and  its  electrical  properties.  A 
second  extension  of  this  technique  would  be  to  exploit  its  potential  multidimensional  search 
capabilities  to  determine  the  frequency  dispersion  of  the  electrical  properties  of  the  skin  and  the 
underlying  breast  tissues.  Finally,  a  long-term  goal  of  this  work  would  be  extension  of  the 
verified  2-D  algorithm  to  three  dimensions. 
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Appendix  A 

Manuscript  submitted  to  IEEE  Antennas  and  Propagation,  August  2002. 

The  manuscript  is  currently  under  revision,  introducing  minor  changes  to  comply  with 
the  reviewers’  comments  prior  to  resubmission. 
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Abstract 

The  work  presented  in  this  paper  is  motivated  by  the  need  to  calibrate  a  new  pulsed- 
microwave  breast  tumor  detection  system  for  patient-specific  skin  parameters.  A  two-dimensional 
time-domain  inverse-scattering  algorithm  based  upon  the  FDTD  method  is  presented  for  deter¬ 
mining  the  skin  thickness  and  the  relative  permittivity  sr-8kin  and  electric  conductivity  askin  in 
the  microwave  range.  The  algorithm  traces  a  search  trajectory  in  the  (er-skin,  <?skin)  parameter 
space.  The  minimal  parameter  estimation  error  along  this  trajectory  yields  a  set  of  approximate 
parameter  values.  It  is  shown  that  the  convergence  of  the  inverse-scattering  technique  depends 
on  the  shape  and  the  duration  of  the  illuminating  electromagnetic  wave  pulse  chosen  for  the 
electrical  parameter  reconstruction.  The  time-domain  nature  of  the  inverse  algorithm  allows  for 
limiting  the  region  of  inversion  using  causality.  Thus,  when  the  parameters  of  the  skin  are  esti¬ 
mated,  the  skin  thickness  can  be  determined  by  comparing  the  measurement  with  a  simulated 
all-skin  response.  Finally,  the  time-domain  inverse-scattering  algorithm  is  tested  for  robustness 
in  the  presence  of  broadband  Gaussian  noise. 
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I.  Introduction 

The  investigations  presented  here  are  motivated  by  the  development  of  a  new  ultra- 
wideband  confocal  microwave  technology  to  detect  and  image  early-stage  breast  cancers 
[1],  [2],  [3],  [4],  [5],  [6],  [7],  [8],  [9],  [10],  [11],  [12]  .  The  new  technology  exploits  the  dielec¬ 
tric  property  contrast  between  normal  breast  tissues  and  malignant  tumors  at  microwave 
frequencies.  Here,  microwave  imaging  is  performed  by  a  planar  antenna  array  contacting 
only  one  side  of  the  breast.  The  antenna  array  elements  collect  impulsive  backscattered 
signals,  which  are  digitally  delayed  to  achieve  coherent  summation  at  the  site  of  a  po¬ 
tential  tumor.  This  scheme  depends  upon  knowledge  of  the  average  dielectric  properties 
of  the  local  breast  tissues.  Patient-specific  calibration  of  the  microwave  imager  requires 
knowledge  of  these  properties. 

In  the  present  research,  a  time-domain  inverse-scattering  technique  is  studied  to  measure 
the  skin  thickness  and  dielectric  parameters  in  the  area  of  the  human  breast.  This  study 
is  motivated  by  several  well-documented  findings  [13],  [14],  [15],  [16].  Skin  thickness  varies 
from  patient  to  patient,  and  also  with  location  on  the  body  of  an  individual  patient.  A 
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number  of  factors  can  cause  thickening  of  mammary  skin.  Although  this  study  is  primarily 
motivated  by  assisting  patient-specific  calibration  of  the  microwave  breast  cancer  imaging 
system,  determining  breast  skin  thickness  can  also  help  to  diagnose  possible  pathologies 
in  the  underlying  tissue  and  in  the  patient  in  general. 

This  paper  first  reviews  the  background  literature  on  dielectric  properties  in  the  mi¬ 
crowave  range,  skin  thickness  in  the  area  of  the  human  breast,  and  numerical  inverse¬ 
scattering  methods.  Next,  calculations  are  presented  to  illustrate  the  importance  of 
knowledge  of  the  correct  skin  thickness  for  determining  the  time  delays  needed  for  the 
image-formation  signal-processing  algorithm.  Then,  results  are  presented  to  show  the 
development  of  a  two-dimensional  time-domain  inverse-scattering  algorithm  for  simulta¬ 
neous  estimation  of  electrical  permittivity  £r-Skin  and  conductivity  crsfcin  of  the  skin  layer. 
This  algorithm  locates  a  search  trajectory  in  the  ( £r-skm ,  Oskin)  parameter  space.  Re¬ 
sults  are  then  presented  which  show  that  the  convergence  of  the  search  trajectory  and  the 
convergence  error  depends  on  the  shape  and  the  duration  of  the  impinging  pulse.  This 
convergence  then  permits  the  skin  thickness  to  be  determined  by  comparing  the  measure¬ 
ment  with  a  simulated  all-skin  response.  Finally,  the  robustness  of  the  algorithm  is  tested 
in  the  presence  of  Gaussian  noise  for  various  signal-to-noise  ratios. 

II.  Background 

A.  Dielectric  Properties  of  Breast  Tissues  in  the  Microwave  Range 

Extensive  dielectric  measurements  up  to  3  GHz  of  both  normal  and  malignant  human 
breast  tissues  have  been  reported  in  the  literature  [17],  [18],  [19],  [20],  [21].  At  6  GHz , 
normal  breast  tissue  dielectric  properties  are  similar  to  those  of  fat  and  vary  in  an  approx¬ 
imate  ±10%  range  about  £r-breast  —  9  for  the  permittivity  and  crbreast  =  0.4  S/m  for  the 
conductivity.  Malignant  tumor  dielectric  properties  are  similar  to  those  of  muscle.  They 
have  nominal  values  of  £T-tumor  =  50  and  atum0r  =  7  S/m.  This  dielectric  contrast  causes 
malignant  tumors  to  have  significantly  larger  microwave  backscattering  cross  sections  rel¬ 
ative  to  normal  tissues  of  comparable  size  and  geometry.  For  breast  skin  tissue,  the  values 
found  in  the  literature  are  approximately  £r-skin  —  36  and  askin  =  4  S/m  in  the  microwave 
frequency  range  [19],  [20],  [21],  [22]. 
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B.  Skin  Properties  and  Thickness  in  the  Area  of  the  Human  Breast 

The  literature  offers  several  studies  on  skin  thickness  in  the  area  of  the  human  breast.  A 
comprehensive  normal  range  of  breast  skin  thickness  and  the  causes  of  thickening  shown  on 
film-screen  mammography  is  given  in  [16].  Breast  skin  thickness  ranges  from  0.7  mm  -  2.7 
mm  depending  upon  the  location.  Specifically,  the  inferior  skin  has  maximum  thickness 
(mean  value  =  1.7  mm);  the  medial  and  superior  skin  thicknesses  both  average  1.5  mm; 
and  the  lateral  skin  has  minimum  thickness  (mean  value  =  1.3  mm). 

The  reported  causes  of  mammary  skin  thickening  are  commonly  categorized  by  their  local 
or  generalized  nature.  The  localized  causes  can  be  carcinoma,  inflammation,  trauma,  fat 
necrosis,  post-biopsy  and  dermatological  conditions.  The  generalized  causes  found  in  the 
literature  are  breast  cancer,  Hodgkin’s  disease,  reticulum  cell  sarcoma,  metastatic  disease, 
radiation  therapy,  inflammation,  surgery,  primary  skin  disorders,  anasarca  and  any  cause 
of  lymphatic  obstruction.  Therefore,  careful  monitoring  of  breast  skin  thickness  can  also 
be  used  as  a  pre-diagnostic  tool  of  the  underlying  tissue  pathologies  or  health  concerns  of 
the  patient  in  general. 

C.  Numerical  Inverse-Scattering  Methods 

Numerical  inverse-scattering  studies  found  in  the  literature  are  based  on  either  frequency- 
or  time-domain  approaches.  With  frequency-domain  algorithms,  the  interaction  of  the 
entire  medium  with  the  incident  field  is  considered  simultaneously  [23],  [24],  [25],  [26].  For 
a  general,  lossy  inhomogeneous  medium  this  can  result  in  a  very  large  number  of  unknown 
parameters.  In  contrast,  time-domain  approaches  can  exploit  causality  to  limit  the  region 
of  inversion,  potentially  reducing  the  number  of  unknowns. 

Time-domain  inverse-scattering  problems  somewhat  related  to  the  present  study  com¬ 
monly  appear  in  the  area  of  geoscience  and  remote  sensing  [27],  [28],  [29].  In  such  prob¬ 
lems,  ground-penetrating  radar  (GPR)  is  used  to  obtain  information  on  subsurface  features 
from  data  collected  over  the  surface.  Signal-processing  algorithms,  including  those  based 
on  FDTD  [27],  [28],  can  be  used  for  target  reconstruction.  Global  inversion  techniques 
have  emerged  involving  the  neural  network  method  [28]  and  the  use  of  genetic  algorithms 
[29].  These  algorithms,  based  on  stochastic  strategies,  offer  advantages  relative  to  local  in- 
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version  algorithms  including  strong  search  ability,  simplicity,  robustness,  and  insensitivity 
to  ill-posedness.  However,  they  require  a  high  number  of  function  evaluations. 

More  relevant  to  the  present  work  is  the  FDTD  /  nonlinear  optimization  techniques  of 
[30]  and  [31].  Reference  [30]  reported  a  one-dimensional  FDTD  formulation  of  an  inverse¬ 
scattering  scheme  for  remote  sensing  of  inhomogeneous  lossy  layered  media.  Here,  a  layer¬ 
stripping  procedure  was  used  to  simultaneously  recover  the  conductivity  and  permittivity 
profiles.  This  work  assumed  an  normally  incident  plane  wave  pulse  and  zero  noise  present 
on  the  received  pulse. 

Reference  [31]  reported  a  two-dimensional  FDTD  inverse-scattering  scheme  for  remote 
sensing  of  the  shape  of  conducting  and  dielectric  targets.  With  this  technique,  as  the 
illuminating  wavefront  sweeps  across  the  target,  causality  is  exploited  to  reconstruct  the 
actual  target  surface  contour  in  a  sequential  and  cumulative  manner.  To  test  the  com¬ 
plexity  of  a  recoverable  structure,  a  plane-wave-illuminated  dielectric  target  with  reentrant 
features  was  reconstructed  from  a  single-point  observation.  Numerical  experiments  tested 
the  performance  of  the  algorithm  in  the  presence  of  noise  and  demonstrated  a  promising 
degree  of  robustness. 

In  the  present  work,  the  one-dimensional  inversion  technique  of  [30]  is  extended  to  two- 
dimensional  illumination  of  a  layered  medium  by  an  infinite  monopole  source  located  at 
the  surface.  Both  the  time  waveform  and  the  duration  of  the  illuminating  pulse  are  studied 
to  explore  optimization  of  the  convergence  of  the  inverse-scattering  technique.  Further, 
the  robustness  of  the  technique  in  the  presence  of  noise  is  examined. 

III.  Single- Parameter  Reconstruction  of  Skin  Dielectric  Properties 
A.  Rationale  and  Motivation 

Fig.  1  shows  a  simplified  2-D  geometry  of  the  17-element  microwave  antenna  array  system 
adjacent  to  the  breast  skin,  underlying  breast  tissue,  and  a  5  —  mm  spherical  tumor  in  the 
breast  tissue  3  cm  below  the  skin  surface  on  the  central  axis  of  the  sensing  antenna  array. 
Signals  received  by  each  antenna  element  of  the  synthetic-aperture  array  are  shifted  in 
time  for  coherent  summation.  We  use  values  of  er-skin  =  36,  £r-breast  =  9  and  principles 
of  geometrical  optics  to  calculate  theoretical  time  delays  between  the  transmitted  signal 
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and  the  backscatter  from  the  tumor  location  for  each  element  of  the  antenna  array  of 
Fig.  1.  Fig.  2  graphs  calculated  time  delays  for  the  extreme  values  of  the  reported  skin 
thickness  range.  We  note  that  the  range  of  variation  (i.e.,  uncertainty)  in  these  delays  is 
comparable  to  the  difference  in  time  delays  required  for  coherent  addition  for  two  adjacent 
antenna  elements  for  each  thickness  value  considered.  This  demonstrates  that  an  incorrect 
assumption  of  skin  thickness  can  lead  to  incoherent  summation  of  the  individual  received 
signals.  Additional  calculations  show  that,  due  to  the  relatively  thin  skin  layer,  the  effect  of 
the  variation  in  £r-Skin  on  the  variation  in  the  calculated  time  delays  is  not  as  pronounced. 
Overall,  the  result  presented  in  Fig.  2  implies  that  we  require  a  patient-specific  calibration 
of  the  confocal  microwave  imager  for  breast  skin  thickness  and  its  dielectric  parameters. 

B.  Basic  Technique 

This  work  focuses  on  the  development  of  a  time-domain  inverse-scattering  algorithm  that 
permits  a  noninvasive  measurement  of  the  near-surface  dielectric  parameters  of  a  lossy 
layered  half-space.  The  new  algorithm  extends  the  iterative  technique  of  [30]  to  the  case 
of  a  2-D  half-space  excited  at  its  surface  by  an  infinitely  long  monopole.  The  principal 
logic  of  this  scheme  is  shown  in  Fig.  3. 

Based  on  an  initial  guess  for  a  set  of  dielectric  parameters  (  sr,  o  ),  the  FDTD  code 
computes  a  trial  backscattered  time  waveform.  The  algorithm  compares  a  trial  signal 
with  the  measured  signal  and  calculates  the  energy-normed  error.  Based  on  this  error,  the 
dielectric  parameters  are  perturbed  by  a  gradient-based  nonlinear  optimization  routine, 
and  the  FDTD  code  computes  a  new  trial  backscattered  signal  based  on  this  new  guess 
for  (er,  a).  This  process  is  repeated  with  varying  eT  and  a  perturbations  as  needed. 

There  are  two  ways  to  exit  the  iterative  loop  of  Fig.  3.  For  exit  1,  the  energy-normed 
error  decreases  below  a  pre-determined  threshold  value.  For  exit  2,  the  error  remains  above 
the  threshold  but  the  estimated  eT  and  a  parameters  oscillate  within  a  narrow  band  of 
values  a  pre-determined  number  of  times.  Both  of  these  exit  criteria  are  chosen  according  to 
physical  meaningfulness  and  to  limitations  imposed  by  available  computational  resources. 
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C.  Results  in  the  Absence  of  Noise 

Here,  we  illustrate  the  convergence  properties  of  the  time-domain  inverse-scattering  algo¬ 
rithm  of  Fig.  3  for  a  skin  half-space  with  only  one  unknown  parameter,  £T-Skin  or  (Tskin.  A 
120-ps  (1/e  full-width)  differentiated  Gaussian  pulse  is  used  as  the  excitation  waveform. 

In  Case  1  (Fig.  4(a)),  we  fix  the  value  of  aSkin  at  the  correct  value  of  4  S/m  while 
iterating  for  £r-skin ■  Here,  the  FDTD  code  computes  the  backscattered  signal  for  a  trial 
value  of  er-skini  compares  it  with  the  reference  backscattered  signal  obtained  for  the  correct 
value  of  er-akini  and  calculates  the  energy-normed  error.  Based  on  this  error,  £r-Skin  is 
perturbed,  and  the  FDTD  code  computes  the  backscattered  signal  for  this  new  guess  for 
the  permittivity.  This  is  repeated  until  one  of  the  criteria  for  exiting  the  iterative  process 
is  met.  Sample  results  are  shown  for  two  initial  guesses  of  £r-skim  99  and  150. 

In  Case  2  (Fig.  4(b)),  we  fix  the  value  of  er-skm  at  the  correct  value  of  36  while  iterating 
for  askin ■  Sample  results  are  shown  for  two  initial  guesses  for  askm,  19  S/m  and  20  S/m. 
In  these  and  similar  simulations,  we  have  found  a  robust  convergence  for  £r-Skin  and  askin 
in  the  absence  of  noise  for  a  variety  of  parameter  perturbation  increments. 

IV.  Two-Parameter  Reconstruction  of  Skin  Dielectric  Properties 
A.  Methods 

This  section  describes  an  optimization  scheme  which  allows  simultaneous  recovery  of  the 
permittivity  and  conductivity  of  the  skin.  This  scheme  generates  a  search  trajectory  in 
the  ( Er-skim  aakin)  space  that  ideally  converges  to  the  reference  values  of  these  parameters. 

The  search  trajectory  is  generated  in  two  possible  ways: 

1.  A  trial  value  of  aSkin  is  assumed  and  a  series  of  FDTD  forward-scattering  runs  is 
conducted  to  find  a  corresponding  value  of  £r-skin  which  minimizes  the  energy-normed 
error  relative  to  the  measured  backscattered  waveform.  This  procedure  is  then  repeated 
for  a  sequence  of  trial  values  of  crSkin  which  vary  uniformly  in  equal  increments  within  a 
±50%  range  about  the  nominal  value  reported  in  the  literature. 

2.  A  trial  value  of  £r-Skin  is  assumed  and  a  series  of  FDTD  forward-scattering  runs  is 
conducted  to  find  a  corresponding  value  of  askin  which  minimizes  the  energy-normed  error 
relative  to  the  measured  backscattered  waveform.  This  procedure  is  then  repeated  for  a 


August  28,  2001 


7 


sequence  of  trial  values  of  £r-skin  which  vary  uniformly  in  equal  increments  within  a  ±50% 
range  about  the  nominal  value  reported  in  the  literature. 

We  here  introduce  the  excitation  waveforms  used  in  our  investigations.  The  results 
section  will  show  that  the  search  trajectory  in  the  (er,  cr)  space  depends  strongly  on 
the  signal  shape  and  duration  of  the  excitation  waveform.  The  three  F^-field  excitation 
waveforms  used  for  the  dielectric  parameter  reconstruction  are  as  follows:  (a)  120  —  ps 
differentiated  Gaussian  pulse;  (b)  10  — ps  differentiated  Gaussian  pulse  with  a  5  —  ps  rise¬ 
time;  and  (c)  5  —  ps  rise  time  ramp  which  has  the  same  maximum  as  the  peak  value  of 
the  10  —  ps  differentiated  Gaussian  pulse. 

It  is  important  to  note  the  following.  The  inverse-scattering  algorithm  processes  only 
a  time  window  TWin(t0W  within  which  it  compares  the  trial  and  measured  signals.  For 
breast-skin  parameter  recovery,  it  is  critical  that  TWindma  be  small  enough  to  avoid  wave 
reflection  from  the  underlying  skin-breast  interface  which  corrupts  the  received  signal. 
However,  Twindow  must  be  large  enough  to  allow  the  measured  backscattered  pulse  to 
contain  meaningful  information  concerning  the  skin  electrical  parameters.  This,  in  turn, 
means  that  only  a  part  of  the  measured  backscattered  waveform  might  be  usable  for 
accurately  generating  a  trajectory  in  the  (£r-skin,  crSkin)  space.  For  example,  in  the  next 
section,  we  will  see  that,  for  the  120  —  ps  differentiated  Gaussian  pulse  excitation  case, 
only  the  leading  edge  of  the  measured  backscattered  pulse  is  used  in  the  reconstruction. 
Conversely,  the  entire  measured  backscattered  pulse  is  used  for  the  10  —  ps  differentiated 
Gaussian  pulse  excitation  case.  Results  will  illustrate  the  significant  consequences  that 
this  has  on  the  calculated  trajectory  in  the  ( £r-akm ,  Vakin)  space. 

B.  Results  in  the  Absence  of  Noise 

B.l  £r—skin  Reconstruction  with  Trial  Values  of  crSkin 

This  section  reports  results  for  the  inverse-scattering  algorithm  iterating  for  er-skin  using 
a  range  of  trial  values  of  aSkin  varying  ±50%  around  the  reference  value  of  askin  =  4  S/m. 
In  all  simulations,  the  time  step  used  in  the  FDTD  code  is  At  =  0.33356  ps. 

Fig.  5(a)  depicts  the  resulting  trajectory  in  the  ( £r-skm ,  &skin)  space  for  the  120  —  ps 
differentiated  Gaussian  pulse  excitation  and  Twind0w  —  400At.  This  trajectory  passes 
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through  the  reference  value  er-skin  =  36.  Fig.  5(a)  also  shows  the  large  error  made  in 
the  er-skin  recovery  if  Twindow  —  600Ai,  which  is  large  enough  to  include  the  unwanted 
reflection  from  the  underlying  skin  -  breast  tissue  interface.  Fig.  5(b)  graphs  the  error 
for  Tendon  for  each  value  in  the  trial  <Jskm  range  for  its  corresponding  estimated  £r-skin  of 
Fig.  5(a).  Fig.  5(b)  shows  that  the  error  function  for  £r-Skin  recovery  using  trial  values 
of  crSkin  suffers  from  local  minima,  and  thus  cannot  be  used  as  a  search  criterion  along 
the  trajectory  in  the  (er-skin,  °skin)  space.  Similar  results  were  observed  for  the  10-ps 
differentiated  Gaussian  pulse  and  5-ps  rise-time  ramp  excitation. 

B.2  a skin  Reconstruction  with  Trial  Values  of  £r-skm 

This  section  reports  results  for  the  inverse-scattering  algorithm  iterating  for  askin  using  a 
range  of  trial  values  of  £r-skin  varying  ±50%  around  the  reference  value  of  £r-skin  =  36 
S/m.  In  all  simulations,  the  time  step  used  in  the  FDTD  code  is  A t  =  0.33356  ps. 

Results  for  the  excitation  signal  waveforms  under  investigation  are  shown  in  Figs.  6,  7 
and  8.  For  the  120  —  ps  differentiated  Gaussian  pulse  (Fig.  6)  and  the  5  —  ps  rise-time 
ramp  (Fig.  7),  Twindow  includes  only  the  leading  edge  of  the  backscattered  pulse  waveform. 
This  yields  in  Fig.  6(a)  and  Fig.  7(a)  linear  search  trajectories  in  the  ( £r-skin ,  &skin)  space. 
For  both  of  these  excitations,  Fig.  6(b)  and  Fig  7(b)  show  a  broad  minimum  of  the  error 
along  the  search  trajectory  as  a  function  of  the  estimated  crSkin  value.  This  minimum  is 
approximately  centered  at  the  location  of  the  assumed  reference  values  of  £r-skin  and  aSk%n ? 
implying  that  minimizing  the  error  could  be  used  as  a  search  criterion  along  the  trajectory 
in  the  ( £r-$kin ,  crs*in)  space.  However,  the  broad  error  minimum  suggests  sensitivity  of  the 
skin  dielectric  parameter  recovery  in  the  presence  of  noise. 

Fig.  8(a)  graphs  the  search  trajectory  in  the  ( £r-skim  askin)  space  obtained  for  the  10— ps 
differentiated  Gaussian  pulse  excitation.  Here,  TWindaw  includes  the  entire  backscattered 
pulse  waveform,  and  the  search  trajectory  shows  a  parabolic  behavior.  We  note  a  sharp 
null  of  the  error  along  the  search  trajectory  at  the  location  of  the  assumed  reference  values 
of  the  skin  dielectric  parameters,  which  suggests  robustness  of  parameter  recovery  in  the 
presence  of  noise.  This  is  explored  next. 
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C.  Results  in  the  Presence  of  Zero-Mean  Gaussian  Noise 

We  now  test  the  robustness  of  the  inverse-scattering  scheme  of  the  previous  section  by 
adding  zero-mean  Gaussian  noise  to  the  simulated  measured  backscattered  signals.  The 
signal-to-noise  ( S/N )  ratio  is  defined  by  20  logw(SignalMAx/NoiseSD),  where  SignalMAx 
is  the  peak  value  of  the  noiseless  backscattered  waveform  and  NoisesD  is  the  standard 
deviation  of  the  zero-mean  Gaussian  noise. 

This  section  presents  results  for  noisy  backscattered  signals  resulting  from  the  10  -  ps 
differentiated  Gaussian  pulse  excitation.  Individual  search  trajectories  in  the  ( er-skin , 
oskin )  space  are  generated  for  five  independent  noise  samples  added  to  the  backscattered 
signal  for  S/N  ratios  of  20  and  40  dB.  Then,  for  each  S/N  level,  the  five  data  sets  are 
averaged  before  input  to  the  inverse-scattering  algorithm  to  yield  a  final  search  trajectory. 
The  latter  procedure  simulates  the  “boxcar  integration”  often  used  as  a  technique  for 
measuring  signals  contaminated  with  zero-mean  noise  wherein  multiple  measurements  are 
averaged  to  enhance  the  deterministic  signal.  For  all  results  shown,  Twinciow  —  80  At, 
where  At  =  0.33356  ps.  Figs.  9(a)  and  10(a)  show  the  generated  search  trajectories 
in  the  (er-skin,  &skin)  space  for  each  noisy  and  the  averaged  backscattered  waveform  for 
S/N  =  20  and  40  dB,  respectively.  As  expected,  for  the  higher  S/N  ratio,  we  observe 
a  smaller  deviation  of  the  individual  trajectories  from  the  trajectory  obtained  using  the 
averaged  backscattered  signal. 

More  importantly,  however,  Figs.  9  and  10  show  that  the  corresponding  error  along  the 
search  trajectory  as  a  function  of  the  estimated  askin  exhibits  a  sharp  null  at  the  correct 
value  of  Er-skin-,  36,  even  for  the  lower  value  of  S/N.  This  unambiguous  indication  of  the 
correct  value  of  £r-skm  despite  a  noisy  background  is  crucial  in  the  determination  of  the 
skin  thickness,  to  be  discussed  next. 

V.  Determining  Skin  Thickness 

We  now  briefly  describe  how  we  can  estimate  the  breast  skin  thickness  once  the  patient- 
specific  value  for  £r-Skin  is  determined.  Note  that  this  procedure  is  relatively  insensitive 
to  crskin.  Here,  the  excitation  signal  shape  and  duration  is  much  less  critical  than  before. 

Due  to  the  dielectric  contrast  between  the  skin  and  the  underlying  breast  tissue,  the 
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propagating  pulse  reflects  off  the  skin  -  breast  tissue  interface.  By  comparing  the  mea¬ 
sured  response  for  the  finite  skin-thickness  case  with  an  FDTD  simulation  which  assumes 
a  homogeneous  skin  half-space  characterized  by  the  £r-$kin  determined  previously,  the  ob¬ 
served  time  of  the  first  reflection  yields  the  data  required  to  determine  the  skin  thickness. 

To  illustrate  this  strategy,  Fig.  11  graphs  the  results  of  four  test  cases  for  an  assumed 
60— ps  differentiated  Gaussian  pulse  excitation  wherein  measured  backscattered  waveforms 
for  skin  thicknesses  of  0.6,  1.2,  1.8  and  2.4  mm  and  £r-skin  —  36  are  simulated  using  2-D 
FDTD  models.  These  waveforms  are  subtracted  from  the  FDTD-calculated  backscattered 
response  for  the  corresponding  skin  half-space.  Then,  the  estimated  time  delay  of  the  peak 
value  of  the  difference  signal  with  respect  to  the  peak  of  the  skin  half-space  response  is 
used  to  estimate  the  skin  thickness.  The  resulting  breast  skin  thickness  values  obtained 
in  this  manner  are  0.51,1.25,1.85  and  2.45  mm,  respectively. 

VI.  Discussion 

A.  Impact  of  the  Excitation  Signal  Shape  and  Duration  Upon  the  Robustness  of  the  Inverse- 
Scattering  Algorithm 

The  results  of  the  previous  sections  show  that  the  robustness  of  the  time-domain  inverse¬ 
scattering  algorithm  for  recovery  of  £r-skin  and  <7skin  depends  on  the  duration  and  shape  of 
the  excitation  waveform.  Specifically,  of  the  excitations  considered,  the  10  —  ps  differenti¬ 
ated  Gaussian  pulse  provides  superior  robustness,  even  relative  to  the  ramp  signal  having 
the  same  rise-time,  if  Twin<iow  includes  the  entire  backscattered  waveform  due  to  the  dif¬ 
ferentiated  Gaussian  pulse.  Apparently,  a  bipolar  excitation  delivers  more  information  to 
the  proposed  inverse-scattering  algorithm  about  the  signal-path  propagation  medium  than 
a  comparably  fast  unipolar  excitation.  The  sharp  error  null  consistently  obtained  along 
the  search  generated  trajectory  in  the  ( £r-akm >  ^skin)  space  for  the  10  —  ps  differentiated 
Gaussian  pulse  excitation,  even  in  the  presence  of  significant  levels  of  zero-mean  Gaussian 
noise,  clearly  illustrates  this  point.  Signal-to-noise  ratios  as  low  as  20  dB  are  adequate  for 
reliable  recovery  of  the  skin  permittivity  using  this  excitation. 
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B.  Linearity  of  Estimated  (rskm  vs.  er-skm  in  the  Case  of  Time-Linear  Excitations 

In  the  previous  sections,  we  observed  that  the  search  trajectory  in  the  ( £T-Skini  ^skin)  space 
is  a  straight  line  for  cases  where  the  leading  edge  of  the  120  —  ps  differentiated  Gaussian 
pulse  or  the  5  —  ps  rise-time  ramp  excitation  was  used  for  simulation  of  backscattered 
signals  and  parameter  reconstruction.  This  was  observed  even  in  the  presence  of  zero- 
mean  Gaussian  noise.  In  this  section,  we  analyze  the  connection  between  the  linear  time- 
dependence  of  the  excitation  pulse  and  the  linearity  of  the  generated  search  trajectory. 

We  begin  the  analysis  with  Ampere’s  law  in  three  dimensions  for  linear,  isotropic,  nondis- 
persive,  lossy  materials: 

^  ^  dE 

V  X.  H  —  ( Jsource  “I"  aE)  +  £  •  (1) 

Here,  H  is  the  magnetic  field  (A/m);  E  is  the  electric  field  ( V/m );  JSOurce  is  an  independent 
source  current  density  (A/m2);  a  is  the  conductivity  (S/m);  and  £  =  £q£t  is  the  electrical 
permittivity  (F/m),  the  product  of  relative  permittivity  £r  with  the  free-space  permittivity 
£q  =  8.854  x  1(T12  F/m. 

In  two  dimensions  and  in  the  absence  of  independent  current  sources,  Equation  1  can 
be  reduced  to  the  TM-mode  scalar  equation  involving  Cartesian  components  of  the  field 


vectors: 


dHy  dH,  dE, 

-m~~m=aE’+e'w 


This  equation  governs  both  the  measured  backscattered  signal  and  the  FDTD-computed 
trial  backscattered  signal: 


dHy>M 

dt 


dHx  ,M  _  TJ!  .  dEz 

-jp-t-.V  +  sr  — 

SH,,r  _  9E,tT 

'^aT-aTE‘'T  +  eT'~aT 


where  the  subscript  M  refers  to  measured  signal  variables  and  the  subscript  T  refers  to 
trial  signal  variables.  The  excitation  for  both  the  measured  and  the  trial  pulse  is  the  same 


and  a  linear  function  of  time: 


Ez.m  —  Ezt  —  at  A  b 


dEzAi 

dt 


dE^r 

dt 
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where  a  and  b  are  constants  with  appropriate  units. 

The  proposed  inverse-scattering  iterative  scheme  computes  the  error  based  on  energy 

difference.  In  essence,  the  iterative  process  attempts  to  minimize  the  following  difference: 

IY  dHyfM  dHXtM  /  dHy,T  dHX£  ..2  _ v  q  /y\ 

u  dt  dt  }  {  dt  dt  n  K  ’ 

In  order  for  this  difference  to  reach  zero,  while  observing  Equations  2,  4,  5  and  6,  the 

following  condition  must  be  fulfilled: 

[(  o’M  ■  (oi  +  b)  +  Em  •  o  )  -  (or  •  +  b)  +  Et  •  u)]2  =  0  (8) 

Clearly,  Equation  8  is  satisfied  when  the  trial  values  of  parameters  are  equal  to  the 

measured  parameter  values,  i.e.,  when  Et  =  £m  and  Ot  =  <7m ■  However,  Equation  8  is 
satisfied  also  for  the  following  condition: 

(at  +  b)  r  (at  +  |  l  ('n'v 

a  ci 

Therefore,  the  minimal  error  condition  of  Equation  8  and  the  final  Equation  9  tells  us  that 
linear  relation  between  et  and  ot  with  appropriate  constants  is  a  sufficient  condition  for 
reconstruction  of  s  and  a  parameters  by  the  proposed  inverse-scattering  scheme  when  the 
excitation  signal  is  a  linear  function  of  time.  This  was  observed  even  for  the  noise  study 
in  the  case  of  the  the  linear,  5-ps  ramp  signal  excitation.  The  potential  benefit  of  the 
demonstrated  linearity  of  the  generated  search  trajectory  in  the  (£r-skin,  &skin)  space  is  as 
follows:  linear  excitation  can  be  used  for  two  trial  sr-$kin  values,  and  the  rest  of  the  values 
along  the  anticipated  linear  ( er-skin ,  Oskin)  trajectory  can  be  obtained  by  interpolation  or 
extrapolation.  Then,  error  values  for  the  sets  of  (£r-skim  a$kin)  obtained  by  this  means 
can  be  tested  with  the  noise-robust  10  —  ps  differentiated  Gaussian  pulse  excitation.  This 
would  significantly  decrease  the  computational  cost  of  the  dielectric  recovery  procedure. 

VII.  Conclusion 


This  work  presented  investigations  of  a  novel  time-domain  inverse-scattering  technique 
that  has  promise  for  permitting  calibration  of  the  pulsed  confocal  microwave  imager  for 
human  breast  cancers.  The  new  technique  is  based  upon  the  use  of  an  FDTD  forward¬ 
scattering  program  element  embedded  within  a  numerical  feedback  loop  containing  a  non¬ 
linear  optimization  routine.  A  systematic  parameter-perturbation  strategy  is  adopted  to 
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permit  searching  of  a  multidimensional  parameter  space  for  the  electrical  properties  of  the 
breast  skin. 

There  are  two  significant  findings  of  this  work.  First,  the  use  of  a  time-linear  (ramp) 
excitation  yields  a  linear  search  trajectory  in  (er-skin,  &skin)  space.  Second,  the  use  of  a 
short  bipolar  excitation  signal  provides  a  sharp  null  of  the  energy-normed  error  along  the 
search  trajectory  even  when  the  backscattered  signal  is  contaminated  with  a  significant 
level  of  additive  zero-mean  Gaussian  noise.  Proper  combination  of  time-linear  and  short 
bipolar  pulse  excitation  data  can  yield  an  efficient  and  robust  search  strategy. 

A  logical  near-term  extension  of  this  work  is  to  construct  an  additional  software-element 
which  would  use  the  same  technique  to  obtain  the  electrical  properties  of  the  breast  tissue 
underlying  the  skin  after  first  determining  the  skin  thickness,  permittivity  and  conductiv¬ 
ity.  A  second  extension  of  this  technique  would  be  to  exploit  its  potential  multidimensional 
search  capabilities  to  determine  the  frequency  dispersion  of  the  electrical  properties  of  the 
skin  and  the  underlying  breast  tissues. 
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Tables  and  Figures 

Fig.  1.  Geometry  showing  the  antenna  array  adjacent  to  the  skin  and  the  path  used  for 
estimating  propagation  time  delays  for  a  5  —  mm  tumor  embedded  in  the  breast  tissue  3 
cm  from  the  skin  surface.  Both  breast  and  skin  tissue  axe  considered  to  be  homogeneous 
with  €f —breast  =  ^  and  Sr— skin  —  36. 

Fig.  2.  Effect  of  skin  thickness  on  propagation  time  delays  for  each  antenna  element  of 
Fig.  1,  assuming  £r-Skin  —  36.  Reported  range  of  normal  skin  thickness  is  0.5  —  2.7  mm. 

Fig.  3.  Time-domain  inverse-scattering  iterative  algorithm  for  recovery  of  dielectric  pa¬ 
rameters  sT  and  a  from  a  measured  (in  this  study,  simulated  measured)  signal.  The  scheme 
exploits  causality  to  limit  the  region  of  inversion. 

Fig.  4.  Inverse-scattering  FDTD  computation  for  a  skin  half-space:  (a)  Convergence  of 
Sr-skin  to  its  correct  value  of  36  for  two  different  initial  guesses  of  sT-skin >  with  oSkin  fixed 
at  the  correct  value  of  4  S/m  .  (b)  Convergence  of  oSkin  to  its  correct  value  of  4  S/m  for 
two  different  initial  guesses  of  aSkin,  with  £r-$kin  fixed  at  the  correct  value  of  36  . 

Fig.  5.  Estimate  of  er-sfem  for  trial  values  of  oskm  for  the  120  —  p.s  differentiated  Gaussian 
pulse  excitation  case,  (a)  Search  trajectory  in  ( £r-skmi  &skin)  space  for  Twind0Vj  =  400At 
and  TWindow  =  600At,  where  A t  =  0.33356  ps.  (b)  Error  vs.  trial  values  of  oskin  for 
Tw indow  ~  400 At. 

Fig.  6.  Estimate  of  oskin  for  trial  values  of  £r-skin  f°r  the  120  —  ps  differentiated  Gaussian 
pulse  excitation,  (a)  Search  trajectory  in  ( sr-skin ,  &skin)  space  for  Twindow  =  400At,  where 
At  =  0.33356  ps.  (b)  Error  vs.  estimated  a$kin  values,  showing  a  broad  null  at  the  location 
of  the  assumed  reference  values  of  the  skin  parameters. 

Fig.  7.  Estimate  of  a$kin  for  trial  values  of  £r-skm  f°r  the  5  —  ps  rise- time  ramp  excitation, 
(a)  Search  trajectory  in  ( eT-Skin ,  ^skin)  space  for  Twindm„  =  lOOAt,  where  At  =  0.33356  ps. 
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(b)  Error  vs.  estimated  askin  values,  showing  a  broad  null  at  the  location  of  the  assumed 
reference  values  of  the  skin  parameters. 

Fig.  8.  Estimate  of  aSkin  for  trial  values  of  £r-skin  for  the  10  —  ps  differentiated  Gaussian 
pulse  excitation,  (a)  Search  trajectory  in  ( £r-skin ,  °skin)  space  for  Twindow  =  80Af,  where 
At  =  0.33356  ps.  (b)  Error  vs.  estimated  oskin  values,  showing  a  sharp  null  at  the  location 
of  the  assumed  reference  values  of  the  skin  parameters. 

Fig.  9.  (a)  Sample  trajectories  in  the  ( eT-skmi  v$kin)  space  for  the  10  —  ps  differentiated 
Gaussian  pulse  excitation  with  S/N  —  20  dB.  (b)  Error  vs.  estimated  oSkm  values  for  the 
averaged-backscattered-waveform  case. 

Fig.  10.  (a)  Sample  trajectories  in  the  ( £r-skin >  &$kin)  space  for  the  10  —  ps  differentiated 
Gaussian  pulse  excitation  with  S/N  =  40  dB.  (b)  Error  vs.  estimated  o3kin  values  for  the 
averaged-backscattered- waveform  case. 

Fig.  11.  Waveforms  illustrating  a  strategy  for  estimating  the  skin  thickness  making  use  of 
the  previously  determined  values  of  £T-skin  and  crSkm-  We  compare  the  simulated  measured 
response  (magnetic  field  component  at  the  location  of  the  antenna)  for  four  finite  skin¬ 
thickness  cases  with  a  simulation  which  assumes  a  homogeneous  skin  half-space. 
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Conference  abstract:  for  the  Bioelectromagnetics  Society  24th  Annual  Meeting,  June  23-27,  2002,  Quebec 
City,  Canada. 


2-D  FDTD  INVERSE-SCATTERING  SCHEME  FOR  DETERMINING  MICROWAVE 
BREAST  SKIN  PROPERTIES:  THE  RESULTS  AND  THE  IMPLIED  OPTIMIZATION 
OPTIONS.  M.  Popovic1  ,  A.  Taflove2*.  1  Department  of  Electrical  and  Computer  Engineering, 
McGill  University,  Montreal,  Quebec,  Canada  H3A  2A7. 2  Electrical  and  Computer  Engineering 
Department,  Northwestern  University,  Evanston,  Illinois  60208,  USA. 

Introduction:  The  work  reported  in  this  paper  is  motivated  by  the  development  of  a  new 
confocal  microwave  technology  to  detect  and  image  early-stage  breast  cancers.  This 
methodology  depends  upon  knowledge  of  the  average  dielectric  properties  of  the  local  breast 
tissues.  Patient-specific  calibration  of  the  microwave  imager  requires  knowledge  of  these 
properties.  Present  results  focus  on  dielectric  properties  and  thickness  of  skin  in  the  area  of 
human  breast.  Well-documented  findings  have  shown  that  skin  thickness  varies  from  patient  to 
patient  and  also  with  location  on  the  body  of  an  individual  patient. 

Objective:  A  three-dimensional  (3-D)  inverse-scattering  finite-difference  time-domain  (FDTD) 
algorithm  would  permit  non-invasive  measurement  of  near-surface  breast  tissue  dielectric 
properties,  based  on  the  information  contained  in  the  backscattered  signal  recorded  with  a  trans¬ 
receiving  antenna  located  on  the  skin  surface.  Initial  investigations  are  done  in  two  dimensions. 
Development  of  a  two-dimensional  (2-D)  FDTD  algorithm  allows  for  study  of  the  inverse¬ 
scattering  scheme  and  the  possibilities  for  its  optimization  with  a  computational  cost  lower  than 
that  of  the  3-D  numerical  investigations. 

Methods:  A  2-D  inverse-scattering  algorithm  based  upon  the  FDTD  method  is  presented  for 
determining  skin  thickness  and  the  relative  permittivity  zr-skm  and  electric  conductivity  askin  in  the 
microwave  range.  This  scheme  applies  an  iterative  technique  to  the  case  of  a  2-D  half-space 
excited  at  its  surface  by  an  infinitely  long  monopole.  The  algorithm  traces  a  search  trajectory  in 
the  (e r-skin,  Qskin)  parameter  space  in  the  following  manner.  A  trial  value  of  Gskin  is  assumed  and  a 
series  of  FDTD  forward-scattering  mns  is  conducted  to  find  a  corresponding  value  of  8 r-skin 
which  minimizes  the  energy-normed  error  relative  to  the  measured  backscattered  signal.  This 
procedure  is  then  repeated  for  a  sequence  of  trial  values  of  <3skm  which  vary  uniformly  in  equal 
increments  within  a  ±50%  range  about  the  nominal  value  reported  in  the  literature.  The  three 
electric-field  excitation  waveforms  used  for  the  dielectric  parameter  reconstruction  are  as 
follows:  (a)  120-ps  differentiated  Gaussian  pulse;  (b)  10-ps  differentiated  Gaussian  pulse  with  a 
5-ps  rise-time;  and  (c)  5-ps  rise-time  ramp  which  has  the  same  maximum  as  the  peak  value  of  the 
10-ps  differentiated  Gaussian  pulse.  The  robustness  of  the  inverse-scattering  scheme  was  tested 
by  adding  zero-mean  Gaussian  noise  to  the  simulated  measured  backscattered  signals  for  signal- 
to-noise  ratios  of  20dB  and  40dB. 

Results  and  Conclusions:  There  are  two  significant  findings  of  this  work.  (1)  The  use  of  a 
time-linear  (ramp)  excitation  yields  a  linear  search  trajectory  in  the  (er.^,n,  <3suf)  parameter  space. 
(2)  The  use  of  a  short  bipolar  excitation  signal  provides  a  sharp  null  of  the  energy-normed  error 
along  the  search  trajectory  even  when  the  backscattered  signal  is  contaminated  with  a  significant 
level  of  additive  zero-mean  Gaussian  noise.  Proper  combination  of  time-linear  and  short  bipolar 
pulse  excitation  data  can  yield  an  efficient  and  robust  search  strategy. 

This  work  was  sponsored  by  the  United  States  Department  of  Defense  pre-doctoral  grant,  under 
Army  Award  number  DAMD17-99-1-9335.  The  authors  would  like  to  thank  Cray  Research  for 
their  computational  resources. 
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Noninvasive  Measurement  of  Breast  Skin  Properties  in  the  Microwave  Range: 
Efficiency  Improvement  of  the  2-D  FDTD  Inverse-Scattering  Scheme 

Prof.  Milica  Popovic 

Department  of  Electrical  and  Computer  Engineering 
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Abstract  -  We  investigate  a  two-dimensional  FDTD  inverse-scattering 
technique  for  noninvasive  determination  of  breast  skin  properties. 
Convergence  of  the  algorithm  depends  on  the  shape  and  the  duration  of 
the  excitation  waveform.  Efficiency  and  robustness  of  this  scheme  are 
improved  by  a  proper  combination  of  time-linear  and  short  bipolar 
pulse  excitation. 

I.  Introduction  and  Background 

Our  work  is  motivated  by  the  need  for  patient-specific 
calibration  of  recently  proposed  systems  [1]  for  microwave 
breast  cancer  detection.  To  obtain  the  tumor  image  within  the 
breast,  the  signals  received  by  a  trans-receiving  antenna  array 
contacting  the  breast  are  summed  coherently  with  a  digital¬ 
processing  scheme  which  assumes  knowledge  of  the  average 
dielectric  properties  of  the  local  breast  tissues.  These 
properties  (relative  permittivity  er_skin ,  conductivity  askin  and 
skin  thickness)  vary  from  patient  to  patient  as  well  as  with 
location  within  the  breast  area  within  an  individual. 

II.  Methods 

Numerical  method .  In  order  to  develop  an  inverse¬ 
scattering  algorithm  which  would  permit  noninvasive 
measurement  of  breast  skin  properties,  we  first  investigate  a 
two-dimensional  (2-D)  time-domain  technique  based  upon  an 
FDTD  forward-scattering  program  element  embedded  within 
a  numerical  feedback  error-minimization  loop.  The  one¬ 
dimensional  inversion  technique  of  [2]  is  extended  to  2-D 
illumination  of  a  layered  medium  by  an  infinite  monopole 
source  located  at  the  surface. 

Simultaneous  recovery  of  the  parameters .  The  scheme 
generates  a  search  trajectory  in  the  (er-Skin ,  CFskin)  space  in  the 
following  manner.  A  trial  value  of  er_skin  is  assumed.  A  series 
of  FDTD  forward-scattering  runs  is  conducted  to  find  a 
corresponding  value  of  askin  which  minimizes  the  energy- 
normed  error  relative  to  the  measured  backscattered  signal. 
This  procedure  is  repeated  for  a  sequence  of  trial  er_skjn  values 
which  vary  uniformly  in  equal  increments  within  a  ±50% 
range  about  the  nominal  value  reported  in  the  literature. 

Excitation  waveforms.  The  time  waveform  and  the 
duration  of  the  illuminating  pulse  are  studied  to  explore 
convergence  optimization  of  the  inverse-scattering  technique. 


Among  the  Ez  -field  excitation  waveforms  under  investigation 
were:  (a)  10-ps  differentiated  Gaussian  pulse  with  a  5-ps  rise 
time;  and  (b)  5-ps  rise  time  ramp. 

III.  Results  and  Discussion 

It  is  shown  that  the  convergence  of  the  inverse-scattering 
technique  depends  on  the  shape  and  the  duration  of  the 
illuminating  electromagnetic  wave  pulse  chosen  for  the 
electric  parameter  reconstruction.  For  the  5-ps  rise  time  ramp, 
the  search  trajectory  in  the  (sr.skin ,  askin)  space  exhibits  linear 
behavior;  the  error  along  the  search  trajectory  has  a  broad 
null.  In  the  case  of  the  10-ps  differentiated  Gaussian  pulse 
(where  the  time  window  of  inversion  includes  the  entire 
bipolar  backscattered  waveform)  the  search  trajectory  shows 
parabolic  behavior;  the  error  along  the  search  trajectory  has  a 
sharp  null  at  the  location  of  the  assumed  reference  values  of 
the  skin  dielectric  parameters,  which  suggests  robustness  of 
parameter  recovery  in  the  presence  of  noise. 

For  the  5-ps  ramp  excitation  waveform,  it  can  be  shown 
that  the  linearity  of  the  search  trajectory  in  the  (er-skin  >  oskin) 
space  follows  from  the  time  linearity  of  the  excitation  signal 
and  the  condition  of  the  minimal  error  condition.  The 
potential  benefit  of  this  linearity  is  as  follows:  linear 
excitation  can  be  used  for  two  trial  er_skin  values,  and  the  rest 
of  the  values  along  the  anticipated  (er_skin  >  crsk;n)  trajectory  can 
be  obtained  by  interpolation  or  extrapolation.  Then,  error 
values  for  the  sets  of  (er_skjn ,  crskjn)  obtained  by  this  means  can 
be  tested  with  the  noise-robust  10-ps  differentiated  Gaussian 
pulse  excitation.  This  would  significantly  decrease  the 
computational  cost  of  the  dielectric  recovery  procedure. 

This  work  was  sponsored  by  the  US  DOD  pre-doctoral  grant 
(Award  Number  DAMD17-99-1 -9335). 
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Dear  Dr.  Popovic: 

On  behalf  of  the  Board  of  Governors  of  McGill  University,  I  write  to  inform  you  that  you  are  appointed  full-time 
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the  University  and  no  modification  thereto  has  any  validity  unless  approved  by  the  Board  of  Governors. 
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University. 
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